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ABSTRACT 

We present a model that investigates the spontaneous emer- 
gence of randomness in equity market microstructure. The 
phase space analysis of our model exposes an endogenous 
source of fluctuation in price and volume. We formulate a 
control problem for maximizing price regularity and stabil- 
ity while minimizing entanglement with the market, rep- 
resenting the NYSE specialists' affirmative obligation to 
maintain 'fair and orderly markets'. 
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1 Overview 

Over the past ten years there has been an unprecedented 
increase in market-making activity, ranging from the glob- 
alized equity markets to progressively more customized 
structures that monetize novel risk dimensions. Tradi- 
tionally the domain of small, specialized, private partner- 
ships, equity market-making in particular has attracted the 
attention of major financial intermediaries, and with it, 
the scrutiny of federal regulators and increasingly risk- 
conscious investors. 

At the same time, a wealth of transaction-level data 
became widely accessible, enabling for the first time ever 
the pursuit of quantitative models for the short-term (high- 
frequency) dynamics of stock prices Q El |3] |4] |5l . The 
interest in such models was further boosted by the apparent 
increase in the frequency and severity of well-publicized 
surprises that have rocked the global financial system re- 
cently, from the collapse of currencies to that of financial 
institutions and the bursting of equity and commodity bub- 
bles. 

Our interest is in the investigation of persistent dynamic 
effects that the rules of market making have on the prices 
of financial assets. Also, we seek to construct a risk man- 
agement framework applicable to equity market-making in 
general, and the monopolist NYSE specialists in particu- 
lar. Our contention is that the microstructure of the inter- 
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actions between market-makers and investors represents an 
endogenous source of randomness 1 6 1 which percolates the 
financial intermediation network and gives rise to empirical 
mesoscopic phenomena that deviate from the predictions of 
the Efficient Mai-kets Hypothesis (EMH). 

The current paper represents a first step in the direc- 
tion of modeling the short-term dynamics of equity mar- 
ket making. Unlike earlier models that focus explicitly on 
the 'microscopic' interactions between buyers and sellers 
Q, our goal here is to attack the problem of the 'ensem- 
ble' dynamics directly. While we believe a complete theory 
will only emerge by exhibiting such ensemble dynamics as 
an appropriate 'thermodynamic' (i.e. large scale) limit of 
some class of microscopic interaction models fsl, this goal 
appears still out of reach. We see our work presented here 
as an attempt to 'bridge' the gap between the burgeoning 
literature on constructive microscopic models and the in- 
creasingly refined empirical studies of 'stylized facts' about 
the equity markets' deviations from the EMH ||9l- 

We begin by presenting our model as a coupled system 
of nonlinear ordinary differential equations (ODEs). The 
state of our system encodes the intra-day price variation 
of an equity, the innovations in its instantaneous trading 
volume, and two slope parameters that couple dynamically 
the fluctuations of price and volume. We proceed to ana- 
lyze this model using techniques from dynamical systems 
theory. Our main result is that, in general, bifurcations in 
phase space lead to endogenous uncertainty in the simul- 
taneous determination of price and volume. We proceed to 
formulate a decision problem faced by the specialist as they 
fulfill their affirmative obligation to the NYSE to maintain 
price regularity while minimizing their interference with 
the market. We end by presenting our next goals in this 
research program. 



2 Model 

The model we present here was motivated by a series of 
papers by Day and collaborators that investigate plausi- 
ble mechanisms that lead to endogenous fluctuations in 



macroeconomics I10II11IIT21 and financial transaction data 
1131 14] 1141 . We proceed by constructing the following set 
of coupled differential equations for the tick-by-tick price 
process xi and the daily trading volume X2'- 
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where 2:3 and 2:4 are unobservable components of the state, 
and z denotes an exogenous input, which can be used to 
model the specialist's control the price setting process. 

Equation Q provides a first-order response dynamics 
for price momentum, which is justified for short enough 
time windows as the result of a Taylor expansion |7J. 
Specifically, one can think of Q as capturing a classical 
'momentum market' when x^ > {) and a 'contrarian mar- 
ket' when 0:4 < 0. Equation (|2ji captures the first order 
condition around the local supply/demand equilibrium'. In 
this context, 2:3 < generally signifies a supply shock, 
while a;3 > signifies a demand shock. Equation (|3} cap- 
tures the second order condition, with 
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This implies for instance that /3i > means a supply shock 
has larger price effects but smaller volume effects than a 
demand shock. On the other hand equation © specifies 
the convexity of the phase picture for the price process xi . 



3 Master Equation 

Dividing and (0) and integrating we obtain 
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From now on we will assume that the initial conditions for 
( Ill4t are such that l33,x^{Q) = /3i 2:4(0). Then (|5j simplifies 

to 2:3/2:4 — /?i//32, transforming (|3 to X2 ~ ii/Ji/Sj'^. On 
the other hand, we see that 

it A = (32X2 

= /9i 2:12:4 — l3iX4z 
= (32X2X4 — (3iX4z 
= X4 {x4 - Piz) (6) 

where the first equality follows by differentiating © and 
the final one by differentiating @ once more. Focusing on 
the homogeneous equation (2 = 0), integration by parts of 
(|6j leads to an integral of motion: 



22:4X4 — 2:4 (2:4 + 2) — Ci 



(7) 



'Observe that, together with 0, this imphes that X2 = X3X1. 



where Ci is the initial value of the left-hand-side. Thus, the 
operator C : f 1-^ 2// — / ( / + 2 1 represents a (nonlin- 
ear) symmetry of the short-term equity dynamics desribed 
by ( Ill4t . Observe that this one-dimensional second-order 
nonlinear ODE applies to any stock, irrespective of its type. 
The parameters /3i and j32 enter in the computation of 2:1 
and X2 from 2:4 and its derivatives: 



Xl 
X2 



/?! ^X4 
(32 X4 



(8) 
(9) 



Letting u = 2:4 and 1; = m = 2:4, we can transform Q into 
a separable ODE by observing that 
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X4 

^ ^ 2:4 = t)-— , 
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which, when substituted into yields 

dv Ci+v{v + 2) 



du 



4 Phase Picture 
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Integrating dlOl l with respect to u and assuming 

V {u — 1) — vq we obtain 
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if Ci < 1 and u ^ A 
if Ci < 1 and u G A 



i+iv) = \v + 1 + C^'\ , 

e^{v) = \v- 

A= [-1 - ^/l-Cl, -1 + VI -Ci), and C2, C3 and 
C4 are the natural normalization constants. 

From the above analysis we can infer the following 
properties of the resulting phase picture: 

1. lim|„|^ookl =oo- 

2. When Ci > 1, the phase diagram consists of two dis- 
connected symmetric components on either side of the 
u = axis (see FigurelSjl. 





Figure 1. Phase diagram for Ci = 0.5 and uq = 1. 
two horizontal lines correspond to — 1 ± ^/2/2. 
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Figure 3. Phase diagram for Ci = and vq = 1. 
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Figure 2. Phase diagram for Ci = —0.5 and vq — \. The 
two horizontal lines correspond to — 1 ± \/6/2. 



3. When Ci ^ 0, the phase diagram crosses the w = 
axis perpendicularly, i.e. 1^1 = cx3 (see Fig- 
ures mj] and |5}. On the other hand when Ci — 0, 

1 4^ I „ = i-u (see Figurellll. 

4. When Ci G (— oo, 1) \ {0}, the phase diagram has a 
unique pole at — 1 + -\/l — Ci (see FiguresQEandl^. 

5. When Ci < 1, the phase diagram is tangent to the 
u = {] axis at —1 — -\/l — Ci (see Figures mjl Hand 

We will use Figure |2] to exhibit the possible solution 
paths. Specifically, we interpret v as the time derivative of 
u and follow the dynamics: 

1. When M > and u > — 1 + V6/2, u explodes to 
+CX). Similarly, when u < and v < —1 — \/6/2, w 
explodes to — oo. 

2. Let B = (0,-l-\/6/2), C = (-0.4,0), D = 
(0.4,0) and E == (O, -1 + \/6/2). These ai-e bifur- 
cation points, in that they entail choices for the path. 
Specifically, points B and E are sensitive to small per- 
turbations along the v axis. Using (|9j we see that these 




Figure 4. Phase diagram for Ci — 1 and vq — 1. 



perturbations measure the value of xi. On the other 
hand, points C and D are sensitive to small pertur- 
bations along the u axis. Using (|9j we see that these 
perturbations measure the value of X2- 
3. Let F denote the point on the phase diagram in the 
fourth quadrant with u — 0.4. Also, let G de- 
note the point on the phase diagram in the second 
quadrant with u — —0.4. Four cycles are possi- 
ble in this phase diagram, depending on choice made 
at the a;2 -bifurcation points: BCEDB, BCGEDB, 
BCEDFB and BCGEDFB. On the other hand the 
a; 1 -bifurcation points can cause explosions, depending 
on the signs of /3i and (32 ■ 

All the versions of the phase diagram (i.e. for any value 
of Ci) exhibit similar bifurcation properties. The resulting 
sensitivity to small perturbations characterizes the system 
dynamics. Furthermore, the resulting path dependence on 
the value of xi and X2 renders it impossible to associate 
a unique value to one if we know the value of the other. 
The residual randomness that is accumulated through the 
time evolution of this system is what we referred to ear- 
lier as an uncertainty principle. Figure |6l shows that the 
absolute value of the correlation p [xi , X2) for Ci = —0.5 
and randomly chosen starting configurations decreases ex- 





Figure 5. Phase diagrams for Ci — 1,2,5 and vq — 1 



Figure 7. Histogram of A2 from a Monte Carlo simulation 
of the dynamics shown in Figure |2] using 1,000 randomly 
selected initial conditions. 



Figure 6. Histogram of absolute value of the correlation 
between xi and X2 from a Monte Carlo simulation of the 
dynamics shown in Figure|2]using 1,000 randomly selected 
initial conditions. 



ponentially at first, and stays strictly bounded away from 1 
(in fact, lower than 0.8). The correlation distributions for 
different values of Ci are qualitatively similar. 
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Figure 8. Normal probability plot of A2 from a Monte Carlo 
simulation of the dynamics shown in Figure fusing 1,000 
randomly selected initial conditions. 



5 Volume and Price Distributions 

We have used the Monte Carlo simulation for the dynam- 
ics in the case Ci = —0.5, /?i = /32 = 1 to compute 
the distribution of volume and price. The tails of each of 
the resulting distributions was fit to a power law and the 
corresponding exponents were recorded. Specifically, for 
each Monte Carlo run i, we computed the empirical mea- 
sure of X2{t) and xi{t) and computed Ai(i) and A2(i) so 
that Pvi (xj < a) ^ a'^^^'^ for j = 1,2. Figure0shows the 
distribution of the resulting exponents A2 (i) for the left tail 
of the volume distribution. The exponents for the right tail 
of the volume distribution exhibit similar qualitative char- 
acteristics (i.e. normality, but with different mean and vari- 
ance). We see that the distribution of A2 is approximately 
normal. Indeed, Figure |S] shows the corresponding nor- 
mal probability plot, which confirms our suspicion that the 



distribution of volume exponents is approximately normal, 
with mean /i2 = —0.88 and standard deviation CT2 = 0.12. 
On the other hand the distribution of the tail exponents 
for prices has a pronounced negative skew, with a mean 
/ii = —0.56 and standard deviation cti = 0.33. In both 
cases, the absolute values of the exponents are substantially 
less than 2, implying light-tailed distributions, both for the 
volume and the price returns. However, these distributions 
are not directly comparable to empirically determined dis- 
tributions. It should be remembered that the randomness 
we observe here arises entirely endogenously, through the 
dynamic bifurcations in phase space, rather than through 
any exogenous inputs. Figures and |9] are shown as ex- 
amples. Different values of Ci, /3i and P2 lead to different 
power distributions, with similar qualitative characteristics. 




Figure 9. Histogram of Ai from a Monte Carlo simulation 
of the dynamics shown in Figure |2l using 1,000 randomly 
selected initial conditions. 



integrated square increments of xi, J (ii) dt. This con- 
straint is often quantified by the exchange as the maximum 
allowable bid-offer spread, which can vary across stocks 
and from time to time, depending on the observed trading 
patterns. 

The specialist's objective as a market maker under the 
NYSE affirmative obligation is the minimization of the 'en- 
tanglement' with the market, which can be measured as the 
absolute value of the correlation between the specialist's 
actions z and the resulting price series xi (see also the dis- 
cussion in 1 13 4 1). On the other hand the specialist's profit, 
n, is generally a multiple of the trading volume, X2, so long 
as the price does not exhibit a pronounced trend. A trend- 
ing market is detrimental to the specialist's profit because it 
often necessitates an inordinate number of one-sided trades 
against the market. Thus, we can model the rate of change 
of n as 



6 Exogenous Input 

Observe that for any exogenous input z, J u{t)z{t)dt = 
J uv~^z{u)du using the fact that du = vdt, where u{t) = 
u ^ z{u) — z{t). If we now return to and add the 
effect of the exogenous input z from (|6j we obtain in lieu 
of MQ\ the following equation: 



dv _Ci+ v{v + 2) + 2/3i / uv- 
du 



2uv 



(11) 



Differentiating with respect to u we obtain after some alge- 
bra: 



^ „ „ I dv\ dv 

UV -—r + UV \ —- ] V— PlUZ = 



du^ 



du J du 



(12) 



Equation (I12> provides us with a tool to choose a control in- 
put z as a function of u in order to attain a desirable (w, v) 
profile. Of course, if we choose z{u) = 0, after some al- 
gebra we recover (|5). On the other hand, if we choose for 
instance z{u) — ku^^, then ( I12> simplifies to the following 
system of two coupled first order ODEs: 

dw vw{l — uw) + k(3i 

du uv^ 

dv 

du 
which can be integrated numerically. 

We are finally in a position to formulate the control 
problem II15I faced by a monopolist market maker at the 
NYSE. In particular, the specialists operate under an 'af- 
firmative' obligation by the exchange to maintain fair and 
orderly markets. Generally, this requirement is interpreted 
as an obligation to complete any one-sided market and min- 
imizing price variation while avoiding as much as possible 
interference with the autonomous evolution of the market. 

The model we developed in this paper allows us to for- 
malize this decision problem. The price regularity con- 
straint can be conceptualized as an upper bound U to the 



n = ciUx2 — C2X1 

= ciUPiP^^xi - C2X1, 



(13) 



for some positive constants ci , C2 > 0. In summary, we 
propose the following control problem for the monopolist 
market maker: 



subject to 
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UV ——^ + UV 
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du 



Xi 
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dv 
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{±iYdt<U 
C2P2X1 - ciPiUxi + (32L < 0, 



(14) 



where L is a lower bound to n from ( I13> that the specialist 
is willing to accept (loss rates higher than \L\ are unaccept- 
able and trigger restructurings that are outside the scope of 
the current model). 

7 Conclusions and Next Steps 



We presented a dynamical model (I1I4> for the short-term 
behavior of stocks. The formal analysis of the model led to 
the determination of an integral of motion, Ci, which de- 
termines the qualitative behavior of the dynamics, as shown 
inFigures[T]El|3]|4|and|5] 

Investigations of the homogeneous phase picture ob- 
tained in jlOt helped us identify bifurcation points, which 
introduce an endogenous source of uncertainty in the high- 
frequency dynamics of price and volume Ill3l l6l. Using this 
framework, we formulated a control problem for the mo- 
nopolist market maker, that attempts to balance the desire 
for price regularity and profit making with the obligation to 
avoid interference with the market ifT4l . 



The solution of the speciaUst's control problem is a nat- 
ural next step of the analysis presented here. This problem 
can be seen from the perspective of the market maker, as 
presented above, in which case we are seeking an optimal 
'market making strategy', z. Alternatively, this problem 
can be seen from the perspective of a risk manager respon- 
sible for the capital adequacy of the market maker. In this 
case it is the maximum allowable loss rate L that is of pri- 
mary interest. Yet a different perspective on this problem 
is that of the regulators who design rules for the exchange. 
In this view, it is the price regularity measure U that is cen- 
tral, as well as the minimum achievable market interference 
(optimal value of the objective function). 

In a different direction, it is desirable to investigate how 
the endogenous randomness introduced by the nonlinear 
dynamics Q interacts with the exogenous input z, and how 
it scales with time. To begin with we may consider the dy- 
namics presented here as a nonlinear filter, taking as an in- 
put random paths z and producing filtered versions. The 
study of this filtering paradigm entails the investigation of 
( I12> . Further, even when z is not random but the result of 
an optimization as shown in ( I14> . it will still interact with 
the randomness produced by the homogeneous dynamics, 
because z effectively continuously readjusts C\ as we can 
see by comparing JlU and dlOt . Finally one would wish 
to recover in some appropriate 'long time' limit the tail ex- 
ponents for the volume and returns distributions that are 
empirically observed ||9||5]E1. 
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